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Abstract 

Jamming transition in traffic fiow (between free and jammed traffic) for liomogeneous 
car following model has been investigated taking into account fluctuations of characteristic 
acceleration/braking time. These fluctuations are defined by Ornstein-Uhlenbeck process. 
The behaviour of the most probable deviation of headway from its optimal value has been 
studied and phase diagram of the system has been calculated for supercritical and subcrit- 
ical regimes of jam formation. It has been found that for the first regime the fluctuations 
of characteristic acceleration/braking time result in coexistence of free moving and jammed 
traffic, that is typical for the first-order phase transition, and in appearance of two steady 
states for the second mode. These states correspond to non-zero values of headway deviation 
at which the formation of jam and congested traffic are possible. Using phase-plain portraits 
method the kinetics of the system transitions has been analyzed for different domains of the 
phase diagram for both regimes. 

Key words. Headway deviation; Lorentz system; Ornstein-Uhlenbeck process; Phase dia- 
gram; Langevin equation 

PACS. 02.50.-r Probability theory, stochastic processes, and statistics; 64.60.-i General stud- 
ies of phase transitions; 89.65.-s Social systems 

1 Introduction 

Nowadays traffic problems attract considerable attention. For a study of the traffic jams formation 
problem thermodynamic stochastic 0, and some hydrodynamic and kinetic theories |P, Q are 
used. These theories are based on car- following model 0, Maxwell model 0, and cellular 
automaton model 0, P). Within the framework of thermodynamic approach jamming transition 
is represented as nonequilibrium ffist-order phase transition. The mentioned method describes 
the deterministic picture of traffic flow. Stochastic theory which is based on the master equation 
allows to find a stationary density of spatially-temporal distribution of traffic jams. However, 
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the fluctuations influence of the parameters characterizing the system is not finally studied. It is 
known that fluctuations do not only play a role of the trigger for the phase transition but also 
lead to essentially new system behaviour, and are a major cause for system's self-organization 
0, 1^. Therefore, at describing the jamming transition the estimation of fluctuations influence of 
the dynamical car characteristics is rational. This is achieved within the framework of synergetic 
approach by the most natural way. 

One of the simplest schemes that describe the self-organization process is the Lorentz system 
[P, y]. At first time it was offered for the atmospheric phenomena description and then used 
for physics, chemistry, biology, sociology problems and so on. Recently synergetic scheme was 



proposed to describe the jamming transition in traffic flow |TD[. In our work offered in Ref. |TD 
the synergetic approach will be developed taking into account internal fluctuations of characteristic 
acceleration/braking time. This parameter is a car characteristic, it indicates a time necessary for 
car to reach the characteristic velocity and plays a key role at the traffic jam formation. We will 
show that internal fluctuations, which have additive noise meaning, lead to the complication of 
traffic flow behaviour. In this work the stationary regime of such system will be considered and 
influence of the characteristic acceleration/braking time in the most probable headway deviation 
from its optimal value will be studied. As it will be shown this headway deviation characterizes 
a phase transition into the state that has traffic jam meaning. 



2 Basic equations 

On the basis of the car following model for one-line highway it is possible to show that dissipative 
dynamics of the homogeneous traffic flow can be represented within the framework of the Lorentz 
scheme To describe the jamming transition we use the synergetic concept of phase transition, 
which is realized as a result of mutually coordinated behaviour of three freedom degrees: an order 
parameter, a conjugate field, and a control parameter. In traffic flow the roles of these quantities 
are played by the absolute value of headway deviation between vehicular from a safety distance h 

r]=\Ax-h\; (1) 

by the deviation of velocity of the rj variation from optimal value h/to — V 

v = Ax- h/to + V, (2) 

where to is the nominal time lag, V is the actual value of the velocity, and by acceleration/braking 
time r, respectively. Here x is the vehicle coordinate. 

Let us consider the simplest homogeneous system with goal of definition of time dependencies 
ri{t), v{t), T{t). To achieve this we will use a phenomenological approach. In the equations of the 
motion it is assumed that in the autonomous mode the evolution of the quantities t], v and r has 
dissipative character and their relaxation to equilibrium values is described by the Debye equation. 
Besides a Le Chatelier principle has a great importance: since the decrease of acceleration/baking 
time r assists to the formation of stable traffic flow, the headway deviation t] and its velocity 
deviation v should vary so that to prevent the growth of r, and as a consequence to impede a jam 
formation. Also the essentially important role is played by positive feedback of values t] and r on 
V. Namely the availability of this feedback is a reason for self-organization that leads to traffic 
jam formation. 
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The Lorentz system takes into account the mentioned above circumstances by the simplest 
way. Taking under consideration fluctuational addition it is determined as follows 

?7 = -J-+V, (3) 
V = -— + g^Tr], (4) 

67; 



Tq-t 



gr7]V + \{t). (5) 



Here the dot means the time differentiation; t^, and t^. are the appropriate relaxation times; 
and g-r are the positive constants. 

System (3)-® represents a base of the self-consistent description for the car-following model. 
Here the first terms on the right-hand side of the equations describe the relaxation of each quantity 
to an equilibrium value. In equation the last term is a usual addition. It is easy to see from this 
equation that in the stationary state headway deviation is proportional to the velocity deviation. 

The second term on the right-hand side of Eq. (^) describes the positive feedback of headway 
deviation r] and acceleration/braking time r on velocity deviation v. This feedback leads to the 
increase of v and is the reason for traffic jam formation. 

Equation (^) differs from (Q), (|) because relaxation of r occurs not to zero but to finite 
value To which represents the time necessary for automobile to reach a characteristic velocity (the 
car property). Minus sign before the last term on the right-hand side of equation may be 
considered as a demonstration of the Le Chatelier principle. 

Quantity A(t) represents an influence of fluctuations of the characteristic acceleration/braking 
time and is defined as Ornstein-Uhlenbeck process: 

(A(t)) = 0, 

mm) = c(t,t') = -e"^, (6) 

where is the noise intensity, t\ is the correlation time of the process A(t). 

Within the framework of the mentioned parametrization the formation of traffic jams is rep- 
resented as a result of spontaneous headway and velocity deviations if the characteristic acceler- 
ation/braking time exceeds a critical value. It is reflected by the appearance of the minimum of 
the effective potential which corresponds to the stationary value of the headway deviation r/o ITOj - 
Therefore we will be interested in 77 evolution further. 

In the general case system have no analytical solution therefore we should use the 

following approximation: 

This condition implies that in the course of evolution the acceleration/braking time r is coordi- 
nated by variation of the headway and velocity deviations. Owing to this condition in equation 
(^) a small parameter can be eliminated that allows us to assume t^f ^ 0. As a result we derive 
expression for the control parameter in the following form 

T = To - grUVV + trX{t). (8) 

To form simpler system let us reduce initial one to the one-parameter model. For that it is 
necessary to express v and r via rj. Equation for v is determined by differentiating with respect 
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to time equation for the velocity deviation v obtained from Eq.(^. Substituting expressions for 
f,-?), and equation @ into and introducing measure scales t^, tj^yi — {(jvdrtrtr]^ , Vm — 
t'^^'^iOvgrtr)'^^'^ 1 Tc = (fl'i;^^)"^) ciud Qvirtfj for time, headway deviation r], velocity deviation v, 
acceleration/braking time r, and for the noise of the characteristic acceleration/braking time Tq, 
respectively, we get: 

ii + fi{l + a + ri'^) = ri{e-(j)-ri^ + ri\{t). (9) 
Here the denotations are introduced: 

e^^. (10) 

Obtained expression allows us to write the evolution equation in the canonical form of the 
motion equation for the nonlinear stochastic oscillator of the van der Pole generator type 

v + i{v)v = f{Tl) + 9{v)Kt). (11) 

where 

7(?7) = l + o + rf, 

fir]) = r]{e-a)-r,\ (12) 
g{ri) = 7]. 

Note that equation (|TTD takes into account reactive behaviour of the system. As is known the task 
of the statistical physics is solved using distribution function P{f], r], t) that represents probability 
density of the availability of the corresponding values of headway deviation rj and its rate of change 
?7 at a given instant of time t. Since jam in the car flow is defined by headway rj and time t we have 
to consider projection of the distribution function in the half-space (//, t). For that the kinetic 
equation for P{r],ri,t) is reduced to the Fokker-Planck equation with respect to P{7],t) function. 



3 Langevin and Fokker-Planck equations 

For derivation of the evolution equation for the probability density P{ri,ri,t) we use continu- 
ity equation for function p{fi, rj, t) which is connected with P by following equality P(?7, //, t) = 
{p{r],ri,t))x. Here ()a stands for averaging over noise A. Continuity equation is constructed by 
standard manner: 

f d - \ d 

\di ^ J ^' ^ -9iv)K't)-QTp{v, V, i), (13) 

where the operator 

d d d 

L{r]^r]) = --f{r])—r] + fi— + f{r])— (14) 
or] orj or) 



is introduced. Averaging equation (|TB|) and using decomposition technique in cumulants []TT], 
we obtain kinetic equation for P: 

I ^ + L{fi, r]) I P{r], r], t) = A{fj, r], t)P{r], r], t). (15) 
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Here 

Q OO 

A{r^,r^,t)=g{^)-j:C^'^L^'\v,r]l (16) 
'^^ k=0 

and moments of the correlation function are defined as follows: 

1 °° 

C^'\t) = -lT',C{t,t-rx)dTx. (17) 

' 

Operators L^'^^ are defined by recurrent formula 

L(^) = [l(^-\l], L^'^ = 9iv)^- (18) 

Here [A, 13] = AB — 13 A is the denotation for the commutator . 

To pass into the half-space [i], t) the distribution function moments are used 

OO 

Pn{v,t)= J{vrP{v,v,t)dv. (19) 



Namely these moments give us an opportunity to find an effective Fokker-Planck equation. Taking 
into account components of integer orders, instead of (|T5p according to we get 

^^Piv,t) = -^D('\v)Piv,t) + ^D^'\v)Piv,t), (20) 

where 

liv) liv) 

D(^)(,) = 4^ (22) 
TiV) 

are the drift and diffusion coefficients, respectively. 

Equation (^) corresponds to the Langevin equation governing the rj evolution 

7) = D«(r^) + v/2C(0)D(2)(r/) ■ at), (23) 
where ^(t) is a white noise with standard properties 

m) = o, mm) = sit-t'). (24) 

For studying the transitions between behaviour regimes of the system let us use path integral 
formalism. To achieve this aim we write Langevin equation in the form of stochastic differential 
equation 

dri = D^^^ {ri)dt + V2CmW)dw, (25) 
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where dw = C,{t)dt represents a Winner process. This notation allows us to get a new process y{t) 

with transition Jacobean dy/drj = (^\^2C^fiUD^^ . Since we use a white noise then for y(t) the 
stochastic differentiation operator can be written 

Constructing by such way the evolution equation for y{t) process we obtain expression for the 
white noise 

= , - , + i (V2C(W))' (27) 

with probability density P{^{t)) oc exp (^—^ J ^'^{t)dt^. Here and further accent means differentia- 
tion with respect to rj. Taking into account the relationship between distributions ^(r/) = P{^)J, 
where J is the Jacobean of transition from ^ to rj field, according to [[14| , |l^, we get following 
expression 

P(?7,r/,t)(xexpf-^ / Cdt) , (28) 



2 

where Onsager-Machlup function C acts as Lagrangian in Euclidean field theory 
Here the effective potential energy U is given by expression 



U 



V2cmW) 2 V ^ 



(30) 



Thus, the system's kinetics will be characterized by Euler-Lagrange equation with respect to the 
function 



4 Stationary states, phase diagrams and phase portraits 

At first let us consider a stationary states. Assuming ?) = in the Euler-Lagrange equation, we 
obtain 

y (/ + V/) - Ic^'hvg' = 0, (31) 

7' (V/^? - Wgf + 2C^'^g'V'g) + C^'^g' [{SWg - 2)^V^ - ^^V^g] = 0. (32) 

Headway deviation stationary values r^o, that correspond to the extremum of the effective potential 
energy function (|30D, are determined as solutions of equation (pi]). The solution of equation (p^) 
gives the point of plateau appearance in [/ vs dependence. Further using definitions for 7(77), 
/(r/), and g{ri), we analyze the transitions for stochastic system in details. The dependence of the 
order parameter stationary values on noise intensity and characteristic control parameter is given 
by equation 



r]^ + rf(l + 2a-e- 2C(^)) - {e - a - 2C(^))(1 + a) + = 0. (33) 
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which are determined by equation (JT^), take 



In this case the correlation function moments (| 
the following form 

CW=5^ C« = 5V,. (34) 

Then from Eqs.(|5BD and ( |5^ ) for the phase diagram curve, that is the boundary of the existence 
region of disordered phase (r/o = 0) corresponding to the free traffic, we get 



6' 



1 



(35) 



First of all we should note that separation of the low-frequency domain in noise spectrum results 
in decrease of the effective time that is necessary for car to reach a characteristic velocity. 

The dependence of the headway deviation stationary values on the characteristic accelera- 



tion/braking time riQ{e) is the solution of equation (|33D and is shown in Fig. la. According to Fig. la 
there are zero minimum in the dependence of effective potential energy on headway deviation, 
which corresponds to disordered state (free moving traffic), and minimum at nonzero stationary 
value ?7o; which meets the ordered state (traffic jam formation). Dotted line in Fig. la shows that 
these minimums are separated by maximum corresponding to unstable state of the system. Thus, 
the increase of the noise intensity results in coexistence of disordered and ordered states inherent 
in first-order phase transition. 

Phase diagram of the system for different values of the noise correlation time tx is represented 
in Fig. 2. There are three domains in it. Domain 1 corresponds to disordered state of the system, 
i.e. free traffic. Domain 2 is characterized by coexistence of disordered and ordered states, i.e. 
free and jammed traffic can coexist at such parameters. Last domain 3 corresponds to ordered 
state of the system. Here only traffic jam exists (minimum of the effective potential energy) while 
free traffic is unstable (maximum of the effective potential energy). 

As it is obvious from Fig. 2 with increase of the noise correlation time the domain corresponding 
to free traffic decreases while domain of traffic jam increases. 

For the consideration of the system kinetics let us use Euler-Lagrange equation 



(36) 



dC _ d dC _dR 
df] dt drj drj ' 

which is supplemented by dissipative function contribution R = ff /2. Its form is typical for gen- 
erating functional method. Taking into account equations (P^D, (PD|), we find differential equation 
of second-order 



V - 



X 



1 fD(2) 



2 D(2) 



.^2 _^ ^(0)^(2)^ _ 



0. 



(37) 



Equation ( p?] ) can be represented as a system of two first-order differential equations. This re- 
organization simplifies the investigation of our model. Such equations set can be solved using 
phase-plane method which allows us to consider the kinetic behaviour of the system on the basis 
of phase portraits in the (?); r/) plane. 
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Phase-plane portraits of such system are shown in Fig. 3 for each domain of phase diagram 
in Fig. 2b for tx = 0.2. Domain 1 has two singular points: D, and unstable node S (Fig.3a). 
Here singular point D corresponds to free traffic flow (disordered state). This point has complex 
character of stability because phase trajectories at ?) < converge to D and diverge from it at 



?7 > |T6|. Coordinate of singular point 5* is determined by solution of equation (^). Solutions 
of equations (|3T|), (|32D for noise intensity C-^-* = 10 are pictured in Fig. 4a. As is obvious from this 
figure the solution of equation ( p2|) is independent on characteristic acceleration/braking time e, 
i.e. it is shown by horizontal line. Thus, coordinate of singular point S for all domains of phase 
diagram may be defined by intersection of given horizontal line with vertical line for appropriate 
e value. 

Phase portrait for domain 2 of the phase diagram is shown in Fig. 3b. It has four singular 
points: D, unstable nodes and O, and saddle S. Similar to the first case singular point D 
corresponds to free traffic. Unstable node corresponds to unstable state of the system, namely, 
to the maximum of the effective potential energy. Saddle S is determined by the solution of 
equation (|32|). Unstable node O characterizes the traffic jam formation (ordered state). 

In Fig. 3c phase portrait meeting the domain 3 is demonstrated. There are three singular 
points: D, saddle S, and unstable node O. Here point D corresponds to free moving traffic 
also, but now phase trajectories which pass through it have closed form at its vicinity |jT6|. As 
previously solution of equation ( P^f ) is presented by saddle 5*. Unstable node O corresponds to 
traffic jam. 

The above consideration represents supercritical regime of the traffic jam formation correspon- 
ding to second-order phase transition. However we see that for given system due to fluctuations 
of the characteristic acceleration/braking time free moving traffic and traffic jam can coexist that 
inherent in the first-order phase transition (subcritical regime). In addition the phase portraits 
(Fig.3a-c) show that singular point O meeting the traffic jam state is unstable. Thus, we should 
consider subcritical regime of the traffic jam formation which is a true reason of self-organization 
and analogous to the first-order phase transition. 

For this let us assume that the relaxation time of headway deviation is the function trj{ri) 
increasing with rj from initial value t^(l -|- /c)~^, > to the final one t^, and is determined by 
the simplest approximation 

, ^ (38) 



where < rjt < 1, k and 7]t are dispersion constant and scale. Then doing the same operations 
as in the supercritical case we obtain equation (|ll]), but now 7(77) and /(//) are defined by the 
following expressions: 



fiv) 



a 



k 



a 1 



1 - v^a^ 

(1+^2/^2)2 

k 

1 + T]^ ja^ 



r^- (1 + 



(39) 



where a = rit/rj^. Further Fokker-Planck equation, Langevin equation, and Onsager-Machlup 
function are determined by the similar manner. However, insertion of Eqs. (|39|) into Eq.(|3lD gives 
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the following expression for the stationary values of headway deviation: 

2C(i)] -m(^e-ad + 2C^^^) + C^°^ = 0, (40) 

where 

k 



d = 1 + 



1 + rj'^/a 



m = l + a + k — ^. 41 

Equations (^41), (|40D and ( PD define the boundary of the existence domain of disordered phase 
(^0 = 0) in phase diagram 

e= , J\, +a{l + k)-26'r,. (42) 

According to this the separation of the low-frequency region in the noise spectrum leads to decrease 
of the time needed for car to reach a characteristic velocity. 

Dependence 77o(^), that is a solution of equation (|^), is pictured in Fig. lb. This figure show 
that increase of the noise intensity C^^^ results in appearance of two stationary states corresponding 
to minimums of the dependence of effective potential energy on headway deviation 77. Hence, we 
can conclude that two stationary values of headway deviation from optimal value exist. At these 
values traffic jam and congested traffic can be formed [1^. The smaller value corresponds to the 
metastable state (dashed curve) while lager value (solid line) meets the stable ordered state of the 
system. These states are divided by unstable state (dotted line) corresponding to maximum of 
the effective potential energy. 

Phase diagrams of the system is shown in Fig. 5 for different values of the noise correlation time 
Tx. Here domain 1 corresponds to disordered state of the system, in other words, only one minimum 
of the dependence of effective potential energy on headway deviation exists meeting zero value. 
This minimum characterizes the free traffic. In domain 2 disordered and ordered states coexist, 
so that in addition to zero minimum the minimum appears meeting the nonzero value of headway 
deviation. This minimum corresponds to traffic jam formation. Domain 3 meets the most complex 
form of effective potential energy. Here the metastable and ordered states of the system coexist. 
Metastable state appearance is caused by displacement of zero minimum (disordered state) along 
the axis of headway deviation of the dependence of effective potential energy. It implies that 
in this domain traffic jam can appear at two different values of headway deviation from optimal 
value. And at the same time zero value of headway deviation corresponds to maximum of effective 
potential energy, i.e. free traffic becomes unstable. Domain 4 characterizes the ordered state of 
the system, i.e. traffic jam at a given headway deviation value and unstable free traffic. Last 
domain 5 describes metastable system's state and is located near the boundary between domain 1 
and 3. At crossing the boundary from domain 1 to domain 5 transition of the displacement type 
from disordered to metastable state occurs. It corresponds to appearance of only one minimum of 
the effective potential energy and means the opportunity of congested traffic formation at small 
value of headway deviation. Free traffic is unstable here too. 

Let us now analyze kinetics of the system using Euler-Lagrange equation (pB]). Phase-plane 
portraits of the system for four domains, that are shown in phase diagram for tx = 0.2 (Fig.Sb), 
are demonstrated in Fig.6. Domain 1 is characterized by the presence of three singular points: 
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stable focus F and saddle S (Fig. 6a). Here singular point D corresponds to free traffic and 
has the same complex character with singular point D for supercritical regime. Coordinates of 
the stable focus F and saddle S are determined by the solutions of equation ( P^ ) with taking into 
account Eqs.(|39D. Solutions of Eqs.([31|), (|32D for subcritical regime for noise intensity C(°) = 10 
are shown in Fig.4b. As is obvious from the figure the solutions of equation (|32D are independent 
on characteristic acceleration/braking time e. Consequently, we may obtain the coordinates of 
singular points F and S by intersection of corresponding horizontal lines with vertical lines for 
appropriate e value. Thus the coordinates of singular points 5, S", F, F' can be defined for all 
phase diagram domains. 

Phase portrait of the system for domain 2 of phase diagram is shown in Fig. 6b. There are 
four singular points: -D, saddles N and O, and stable focus F' . As in the Fig.6a singular point D 
corresponds to free traffic, but now phase trajectories which pass through it have closed form at its 
vicinity ||16[. Point F' is determined by the solution of equation (p2D. Saddle O characterizes the 



traffic jam formation. Saddle N corresponds to unstable state of the system, namely, to maximum 
of the dependence of effective potential energy on headway deviation rj. 

In Fig. 6c phase portrait for the most complex domain 3 of phase diagram is demonstrated. 
The maximal number of singular points is placed here: saddles M, N and O, stable focus F, 
and stable node F' . As in the Fig.6a coordinates of points F and F' are determined by solutions 
of equation (P^. Singular points and O have the same meaning as in Fig. 6b. Saddle M 
corresponds to metastable state of the system, namely, to the formation of congested traffic. 

Phase portrait in Fig.6d corresponds to domain 4 of phase diagram. As is shown above the 
solutions of equation ( p2D are represented by saddle S" and stable node F' . Saddle O corresponds 
to jammed traffic. 



5 Discussion 

The above consideration shows that for supercritical regime, which corresponds to the second- 
order phase transition, due to fluctuations of the characteristic acceleration/braking time free 
and jammed traffic can coexist. As is known such a picture is inherent in first-order phase 
transition. For subcritical regime, which meets the first-order phase transition, due to regulation 
of fluctuations intensity of characteristic acceleration/braking time a given system can be passed 
from free traffic regime to traffic jam formation. The last state can appear at different values of 
headway deviation between transport facilities. State of the system corresponding to the smaller 
value of headway deviation may be represented as a congested traffic regime. However, the correct 
answer can be defined using additional investigation of dynamical car characteristics, but such 
study is out of our consideration framework. 

It is shown that transition to low-frequency fluctuations of acceleration/braking time (increase 
of Tx) leads to oppose to free moving traffic regime bounding its domain. It is worthwhile to note 
that for subcritical regime of self-organization the singular points, that correspond to effective 
potential energy extremum, are saddles. Such situation is typical for description of stochastic 



system by Euclidean field theory method |17 . 
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Captions 

to the paper "Jamming transition with fluctuations of characteristic acceleration/braking time 
within Lorentz model" by A.V. Khomenko, D.O. Kharchenko, O.V. Yushchenko 

Fig.l - Dependence of the stationary value of headway deviation rjQ on characteristic accelera- 
tion/braking time € = Tq/tc at different noise intensities C^'^^ of e (curve 1 corresponds to C^^^ = 0, 
2 - C(o) = 5, 3 - C(°) = 10) for a = 1 and tx = 0.2: (a) for supercritical regime; (b) for subcritical 
regime at a = 0.1 and k = 1. 

Fig. 2 - Phase diagrams of the system at supercritical regime for different values of noise 
correlation time at a = 1: (a) tx=0.0; (b) tx = 0.2; (c) tx = 0.4. Domain 1 corresponds to free 
moving traffic, 2 - traffic jam and free traffic, 3 - traffic jam. 

Fig. 3 - Phase portraits of the system at supercritical regime for different domains of phase 
diagram at C^^^ = 10, a = 1, and tx = 0.2: (a) domain 1 at e = 0.5; (b) 2 - e = 1.5; (c) 3 - e = 3. 

Fig. 4 - Dependence of solutions rjo of equations ( ^1) and (|32|) for stationary states on charac- 
teristic acceleration/braking time e = tq/tc at C^^^ = 10, tx = 0.2, and a = 1: (a) for supercritical 
regime; (b) for subcritical regime at a = 0.1, A; = 1. 

Fig. 5 - Phase diagrams of the system at subcritical regime for different values of noise cor- 
relation time at a = 0.1, 0" = 1, and k = 1: (a) tx=0.0; (b) tx = 0.2; (c) tx = 0.4. Domain 1 
corresponds to free moving traffic, 2 - traffic jam and free traffic, 3 - metastable congested traffic 
and traffic jam, 4 - traffic jam, and 5 - metastable congested traffic. 

Fig. 6 - Phase portraits of the system at subcritical regime for different domains of phase 
diagram at a = 0.1, a = I, k = 1, and tx = 0.2: (a) domain 1 at C^^-* = 10, e = 1; (b) 2 - 
C(o) = 2, £ = 1.5; (c) 3 - C(°) = 10, e = 2; (d) 4 - C(°) = 10, e = 3. 
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